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Abstract. Three-dimensional hydrodynamic calculations are performed in order to investigate mass transfer in a close binary 
^ I system, in which one component undergoes mass loss through a wind. The mass ratio is assumed to be unity. The radius of the 

mass-losing star is taken to be about a quarter of the separation between the two stars. Calculations are performed for gases 
with a ratio of specific heats 7 = 1.01 and 5/3. Mass loss is assumed to be thermally driven so that the other parameter is the 
sound speed of the gas on the mass-losing star. 

Here, we focus our attention on two features: flow patterns and mass accretion ratio, which we define as the ratio of the mass 
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' accretion rate onto the companion. Mace, to the mass loss rate from the mass-losing primary star, Mioss- We characterize the 

, flow by the mean normal velocity of wind on the critical Roche surface of the mass-losing star, Vr. When Vr < OAAfl, 

where A and Q are the separation between the two stars and the angular orbital frequency of the binary, respectively, we obtain 
Roche-lobe over-flow (RLOF), while for Vr > 0.7AQ, we observe wind accretion. We find very complex flow patterns in 
between these two extreme cases. We derive an empirical formula of the mass accretion ratio as 0.18 x iq-O-TsYr/asi -^^ 
Q ' low velocity regime and 0.05 [VR/AQy^ in the high velocity regime. 

C/3 . Key words, wind accretion - hydrodynarmcs simulation - binary stars 
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1. Introduction are usually much lower than those required to maintain the 

5; symbiotic activity of the systems. For example, Nussbaumer 

^ Wind accretion plays an essential role especially in the evolu- (,1991, pointed out that in most symbiotic stars the typical rates 

- - tion of detached binary systems such as symbiotic stars, pre- ^-^^^ j^ass loss red giant are in the range 1 x IQ-HIq yr^i 

cursor of peculiar red giants, C Auriga stars and massive X-ray ^ ^ IQ-HIq yr"! and that the expected accretion efficiency 

binaries. In the majority of symbiotic stars, the system con- ^j^^ ^^p^^^g probably no better than one percent. Hence, 

tains a hot component - beheved to be a mass-gaining white ^^e accretion rates associated with most symbiotic white dwarfs 

dwarf - and a cool component, the mass-losing red giant (e.g., ^^.g probably in the range of 1 x IQ-HIq yr'^ to 1 x IQ-^^Mq 

Mikolajewskai997 for a review). Recent studies indicate that yj.-i^ ^^ich are too low to power the typical luminosities of 

many of the red giant components in symbiotic stars do not fill symbiotics (> IOOOLq) unless they are all in a late phase of hy- 

their Roche lobe, i.e., Rno ^ ^i/2, where Rno is the radius ^j^ogen shell-flash (novae) on a white dwarf (e.g., Sion & Ready 

of the red giant component and is the distance of the inner ^gg^ gjon & Staii-fieId'T994^. It will nevertheless be shown in 

Lagrangian point from flie centre of the red giant (e.g., Miirset t^e present paper that the mass accretion rate is smaller than 

& Schmid 19991. This possibly discloses that symbiotic stars ^^at expected by a simple Bondi-Hoyle picture, 
are - with probably only the exception of T CrB (e.g., Anupama 

& Mikolajewska! 1999 1 - well detached binary systems. Thus, It has recently been suggested that some of the mass- 
mass transfer in symbiotic stars seems to be driven by wind, accreting white dwarfs in symbiotic stars are the progenitors of 
and wind mass loss is therefore a key ingredient for triggering Type la supernovae (e.g., Hachisu & Kato 2001 1. On the other 
symbiotic activity on a white dwarf companion. On the other hand, it has long been discussed that the efficiency, the ratio of 
hand, the mass accretion rates expected in the Bondi-Hoyle pic- the mass captured by the white dwarf to the mass lost by the 
ture (Bondi & Hoyle .1944; for such detached binary systems red giant, is very small (say, a few per cent) and not enough to 
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increase the white dwarf mass to the Chandrasekhar mass Hmit 
(e.g., Kenyon et al. 1993 ). 

Wind velocities in cool red giants are typically a few times 
10 km s^^ and comparable to the orbital velocity in most sym- 
biotic binary systems. Therefore, flow patterns are probably 
something between Roche lobe overflow (RLOF) and Bondi- 
Hoyle wind accretion kinds of flows. The efficiency of capture 
in symbiotics should therefore not be estimated by the Bondi- 
Hoyle picture but by direct simulations of intermediate, com- 
plicated flow patterns. 

Symbiotic stars are not the only systems in which wind ac- 
cretion play a significant role. The most adopted model for the 
formation of peculiar red giants. Barium, CH and S stars, re- 
quire a binary system containing an Asymptotic Giant Branch 
star (AGB) transferring mass via its stellar wind to a main se- 
quence companion which becomes polluted in Carbon and s- 
process elements (see e.g. Boffin & Jorissen 1988 1. The AGB 
then evolve into a white dwarf while the companion will appear 
as Carbon or s-process rich and, when on the giant branch, as 
a peculiar red giant. Boffin & Zacs ( 1994i have shown that an 
accretion efficiency of only a few percent is enough to explain 
present Barium stars. Objects related to these are the extrinsic 
S stars which also show symbiotic activity (e.g. Carquillat et 
al.|1998i and bipolar planetary nebulae (e.g. Mastrodemos & 
Morris 1998). 

( Aurigae systems are another kind of system where wind 
accretion plays a role. In these eclipsing double-lined binaries, 
a main sequence star has an accretion wake produced by the 
wind of its K supergiant companion. The eponymous system, ( 
Aur, is a binary with a 972 days orbital period and containing a 
5.8 Mq K4 lb and a 4.8 Mq K B5 V star (Bennett et al. l996 l. 

Finally, accretion wakes are also reported for high mass X- 
ray binaries (HMXB). These systems consist of a compact ob- 
ject, neutron star or black hole, orbiting a massive OB primary 
star which has a strong stellar wind. The X-ray emission is be- 
lieved to be due to accretion of matter on the compact compan- 
ion. 

There have been already many numerical studies of stel- 
lar wind in a close binary system. Biermann J1971I I com- 
puted two-dimensional stellar wind using a characteristic 
method. Sorensen, Matsuda & Sakurai J1975> performed two- 
dimensional finite difference calculation of stellar wind emit- 
ted from a Roche lobe filling secondary. They used the Fluid 
in Cell (FLIC) method with first order of accuracy and a 
Cartesian grid. They computed Roche lobe overflow (RLOF) 
as well. Sawada, Hachisu & Matsuda (1984 1 calculated two- 
dimensional stellar wind from a contact binary using a Beam- 
Warming time implicit finite difference scheme and a gener- 
alized curvilinear coordinate. Their main goal was to com- 
pute the angular momentum loss rate from the system, which 
is an important factor to define the evolution of the binary 
system. Sawada, Matsuda & Hachisu (1986) computed two- 
dimensional stellar wind from a semi-detached binary sys- 
tem using the Osher upwind scheme and a generalized curvi- 
linear coordinate. They observed a transition from RLOF to 
stellar wind by increasing the wind speed on the mass-losing 
star. The present study is a three-dimensional version of their 
study. Matsuda, Inoue & Sawada, (.1987j also conducted a 



similar study and discovered a flip-flop instability of the bow 
shock formed around a compact mass accreting object (see also 
Boffin & Anzer, 119941 and on the stability issue, Foglizzo & 
Ruffert, 1997 and 1999, and Pogolerov, Ohsugi & Matsuda, 
2000). However, it was found that the phenomenon was charac- 
teristic to two-dimensional case (Matsuda et al., 1992 Ruffert, 
il996j . In the present work we perform three-dimensional cal- 
culation in which we do not observe the flip-flop instability. 

These simulations concerned the case of a compact object 
moving in a wind and were mostly aimed at estimating the va- 
lidity of the Bondi-Hoyle accretion rate. Binary effects were 
partially simulated by including velocity or density gradient in 
the wind. There are a few simulations however which simulate 
in full mass transfer by wind in a binary system (Blondin et al. 
1991, Theuns & Jorissen 1993, Blondin & Woo 1995, Theuns 
et al. 1996, Walder 1997, Mastrodemos & Morris 1998, Dumm 
et al. 2000, Boroson et al. 2001, Gawryszczak, Mikolajewska, 
& Rozyczka 2002). Theuns & Jorissen (J993) and Theuns, 
Boffin & Jorissen ( 199^ TBJ96 in the following) performed 
three-dimensional hydrodynamic simulations of a 3 Mq AGB 
transferring mass through its stellar wind to a 1 .5 Mq main se- 
quence companion, in order to test the wind accretion model 
for the origin of peculiar red giants. They performed simula- 
tions using a polytropic equation of state with 7 =1., 1.1 and 
1.5 and found mass accretion ratios (that is, mass accretion 
rate/mass loss rate) of 1-2% for the j ~ 1.5 case and 8 % for 
the other models, i.e. about ten times smaller than the theoreti- 
cal Bondi-Hoyle estimate. They also observed that this value is 
dependent on resolution. 

Walder (1997) presented simulations of wind accretion in 
well separated binaries for three different cases: a HMXB, ( 
Aur and a barium star progenitor. He obtains mass accretion 
ratio (mass accretion rate/mass loss rate) of 0.6, 3 and 6 %, re- 
spectively. Walder & Folini (2000) also presented a nice review 
of wind dynamics in symbiotics binaries, with an emphasis on 
the influence of the radiation field of the accretor. 

More recently, Dumm et al. ( I2000> presented three- 
dimensional Eulerian isothermal simulations in order to repre- 
sent the symbiotic system RW Hya which seems to present an 
accretion wake, as indicated by the reduced UV flux observed 
at phase <I> = 0.78. They found that 6% of the M-giant wind is 
captured by the companion. 

In this paper we perform three-dimensional numerical sim- 
ulations of mass transfer by wind in a binary system with a 
mass ratio of unity. The full gravitational forces of the two com- 
ponents are taken into account. This is therefore different from 
the e.g. TBJ96 or Dumm et al. (2000) studies, where in order 
to account for the not very well known acceleration mechanism 
of the wind, the gravitational force of the primary is partially 
or totally reduced. 

Here, we will concentrate on the derivation of the mass ac- 
cretion ratio. We will leave for a following paper the physical 
discussion of the flow, including a discussion of angular mo- 
mentum loss from the system. 
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2. Assumptions and numerical methiod 

2.1. Model 

Consider a detached binary system of equal mass: one larger 
star blows stellar wind which the other, more compact, compo- 
nent partly accretes. We normalize the length by the separation 
of the two stars. A, and the time by 1 /fi, where Q, is the angular 
velocity of the binary system. The surface of the mass-losing 
star is assumed to be an equi-potential surface, whose mean ra- 
dius is 0.25. The radius of the mass accreting star is assumed 
to be 0.015 or smaller 

From an OB star, stellar wind is mainly accelerated by line 
absorption of UV radiation. In symbiotic stars or in precursors 
of peculiar red giants (PPRGs), the stellar wind mechanism 
from the cool star is poorly known, and may imply radiation 
onto dust grains. One needs however to still bring the matter 
far enough away from the star so that dust grains can form. In 
more evolved red giants, this may be done through the effect 
of pulsation. It is however still not clear what the mechanism 
could be in non-pulsating red giants, as they seem to exist in 
some symbiotic or extrinsic S stars. 

It is therefore clear that simplifying assumption have to be 
made. In the present study, we will thus simply assume that the 
wind is generated by thermal pressure, like in the solar wind. 
More complicated wind structures will be studied in a forth- 
coming paper. 

We assume that the gas is an ideal one and is character- 
ized by the ratio of specific heats, 7. In the present work we 
consider two extreme cases : 7 = 5/3 and 7 = 1.01. In the lat- 
ter case, the gas is almost isothermal except where there are 
shocks. This isothermality may be explained as an outcome 
of good thermal conduction. We neglect other complex effects 
like magnetic fields, radiation and viscosity (except numerical 
one). 

We have carried out our simulations for the case of a mass 
ratio unity. Although this may not be representative for PPRGs, 
symbiotic stars or HMXBs, this should not change very much 
our results. 

2.2. Method of calculation 

We solve the three-dimensional Euler equations using a finite 
volume scheme. As Riemann solver, we use the SFS scheme 
and we apply the MUSCL method to interpolate physical vari- 
ables in a cell; the details on the scheme is described in Makita, 
Miyawaki & Matsuda (2000), and RLOF simulations using the 
same method are given in Fujiwara et al. (2001 ). 

Calculations are done in the rotating frame of a binary and 
include the gravity of the two stars, the centrifugal and Coriolis 
forces. We use a Cartesian coordinate with the origin at the 
mass accreting star. The computational region is —1.5 < x < 
0.5, —0.5 < y < 0.5, < z < 0.5 and we assume symmetry 
around the orbital plane (Fig.^. 

The region is divided into 201 x 101 x 51 cells. The mass 
accreting companion star is represented by a cubic hole with 
cells of 3 X 3 X 3, including the lower half of the hole under 



the orbital plane. Its physical dimension is 0.03 x 0.03 x 0.03 
in our units. 

In order to study in details the structure of the flow near the 
mass accreting object, we introduce a three levels nested grid. 
At the second level, the region -0.25 < x < 0.25, -0.25 < 
y < 0.25 and < z < 0.25 is further divided into 102 x 101 x 
51 cefls. At the third level, the region -0.125 < a; < 0.125, 
-0.125 <y < 0.125 and < z < 0.125 is once more divided 
into 104 X 104 X 54 cells. The cell size of the third level is thus 
1/4 of the original first level cell. The hole of 3 x 3 x 3, in the 
first level, is represented by a hole with 12x12x12 cells in 
the third level. 

A typical symbiotic star, ( Auriga or PPRG has an orbital 
period of a few years and a separation of a few AU. In this 
case, the cell size of 0.03 is of the order or even larger than 
a main sequence star. This is thus applicable for C Aurigae or 
PPRGs. For S or symbiotic stars, we should simulate the case 
of an accreting white dwarf. This would require a hole about 
100 times smaller This is well beyond present days numerical 
capabilities. This is also true if one would like to simulate the 
compact accreting object in HMXBs. 

One should note however that the size of our accreting ob- 
ject is still smaller than the one adopted by e.g. TBJ96. They 
have shown that the mass accretion rate decreases with de- 
creasing accreting size and our value is similar to their high- 
est resolution model. We thus also consider a case of a hole 
having a size of 3 x 3 x 3 cells (i.e. of physical dimensions 
~ 0.01 X 0.01 X 0.01 at the third level in order to investigate 
the effect of the size of the mass accreting object on the results. 



z 




Fig. 1. Computational region. 



2.3. Boundary conditions and initial condition 

The outer parts of the mass-losing star is filled by a gas of den- 
sity (0 = 1 and sound speed Cg. If the pressure inside the mass- 
losing star is larger than that outside the mass-losing star, gas is 
ejected from the surface. The velocity of the ejected gas is com- 
puted by solving a Riemann problem between the two states. 

At t = the computational region except the mass-losing 
star is filled by a tenuous gas with density po = 10^^, velocity 
uq ^ vq — Wo — and pressure po = 10^''/7. The com- 
putational region is gradually filled with the gas ejected from 
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the mass-losing star. The outside of the computational region is 
filled by the initial gas all the time. Mass outflow/inflow from 
the outer boundary can occur; the amount of outflow/inflow 
is also calculated by solving Riemann problems between two 
states. This (outside) boundary condition, which we call the 
ambient condition, insures the stability of the computation. The 
inside of the mass accreting hole is almost vacuum and the gas 
approaching the hole is absorbed. 

We start the calculation at f = and terminate it when 
the system becomes steady. In the wind cases half an orbital 
period is enough to reach a steady state, but in the RLOF cases 
we have to calculate at least a few orbital periods. 

3. Numerical results 

3.1. Parameters 

As was described earlier, we adopt 7 and the dimensionless 
sound speed of gas in the mass-losing star, Cg, as the model 
parameters. In the present work, we consider the cases of 7 = 
1.01 and 7 = 5/3. When Cg is less than some critical value, we 
expect RLOF type flows, while in the other cases, we expect 
stellar wind kind of flows. As is shown in Tabled we study the 
cases of Cs = 0.5 — 2.5 for 7 = 1.01, and Cs = 0.75 — 2.5 for 
7 = 5/3. 

In PPRGs and symbiotic stars, the sound speed at the sur- 
face of the mass-losing giant is of the order of 5-20 kms^^ 
while the orbital velocity is of the order of 10-40 kms^^. In ( 
Aurigae, both these values are around 40-100 kms^^, while in 
HMXBs, the wind speed is of the order of a few thousands of 
kms~^ while the orbital velocity is several hundreds of kms~^. 

The parameter Cg has no physical significance, because the 
mean radius, 0.25, of the mass losing star is chosen arbitrar- 
ily for numerical purposes. Using the results of our numerical 
simulations, we therefore calculate the mean value of the ver- 
tical component of the gas speed on the critical Roche surface, 
Vr. We classify the flow pattern using Vr, which is shown in 
Table [2 as well. It has to be emphasized that Vr is not con- 
stant along the critical Roche surface, especially in the Roche 
lobe type of flows. For these, it may be better to use as a 
parameter. For the wind accretion type of flow, which are more 
isotropic, Vr is more appropriate. 



We define the flow to be of the RLOF type if the flow circulates 
around the accreting object in counterclockwise direction. 

If we raise Cg further, the gas escapes from the mass los- 
ing star in all directions, and part of the gas rotates in a clock- 
wise direction around the companion. These two flows, rotat- 
ing counterclockwise and clockwise, collide behind (i.e. to the 
right) the accreting object to form a bow shock. This is a typical 
wind accretion flow. 

In the B-D cases for 7 = 5/3 and the a case for 7 = 1.01, 
we see from Table^that Vr < 0.4, i.e. Vr is smaller than the 
escape velocity. The flow is of the RLOF type. Similar types of 
flows have been obtained by other authors, e.g. Bisikalo et al. 
(1998). 

Figurel^shows the density contours for model B (7 = 5/3, 
Vr = 0.01), and Figure|5]depicts the streamlines for the same 
case. In this case, since 7 = 5/3, the temperature of the disc 
is so hot that the mass accretion ratio /, which is defined in 
Sect, is not 1 but 0.34. Spiral shocks, typical of RLOF 
flows (Makita et al. .2000. Matsuda et al. .2000. Fujiwara et al. 
2001 ) are not seen in this case. 

One should note that contrarily to the impression one might 
have when just looking at the figures, we have checked that 
there is no inflow from the outer boundaries in these simula- 
tions. 

Figure |6l shows the comparison of model a (7 = 1.01, 
Vr = 0.06) and D (7 = 5/3, Vr = 0.12). We observe spi- 
ral shocks, although with a very much distorted shape. 




I -0.5 
.5 



-1.5 ^ 0.5 



-5.04 -3.36 -168 0.00 

Fig. 2. Density contours in the orbital plane of the typical 
RLOF case, but in a semi-detached binary system (7 = 1.01, 
Cs = 0.02). 



3.2. Roche lobe overflow 

In a typical Roche lobe overflow (RLOF), gas on the surface of 
the companion flows through the L 1 point and forms a narrow 
LI stream. The LI stream deflects its motion to the negative y 
direction in the present configuration due to the Colioris force, 
and then circulates counterclockwise about the mass-accreting 
star to form an accretion disc. We show an example of such 
typical ROLF flow in the case of 7 = 1.01 in Fig. Hand |3] 
Here we assumed the star to fill its Roche lobe, so this is not 
the main aim of this paper. They will serve as comparison for 
the B-D and a cases discussed below. 

If the temperature of the surface of the companion star, c^, 
is raised, the LI stream becomes thicker (Fujiwara et al.. .2001j . 




Fig. 3. Streamlines and Mach number counters of the typical 
RLOF case (7 = 1.01, Cs = 0.02). 



Nagae T. et al.: Wind accretion in binary stars 



5 



Table 1. Model parameters. The parameters Cs, Vr and / are the sound speed on the surface of the mass losing star, the normal 
velocity at the critical Roche surface and the mass accretion ratio. 





7 = 5/3 


7 = 1.01 


Nested Grid (7 = 


1.01) 


Cs 


Case Vr f 


Case Vr f 


Case Vr 


/ 



0.5 
0.75 
0.9 
1.0 
1.2 
1.4 
1.6 
1.8 
2.0 
2.5 



B 
C 
D 
E 
F 
G 
H 
I 
J 



0.01 
0.03 
0.12 
0.68 
1.49 
1.86 
2.20 
2.52 
3.29 



3.4 X 10"^ 

1.0 X 10"^ 
1.9 X 10"^ 

3.3 X 10~^ 

1.1 X 10"^ 

7.4 X 10"^ 
3.6 X 10"^ 
2.1 X 10"^ 
7.9 X 10""^ 



0.06 
0.46 
1.09 
1.41 
1.95 
2.42 
2.86 
3.29 
3.71 
4.75 



1.7 X 10"^ 
9.2 X 10"^ 
2.4 X 10"^ 
1.1 X 10"^ 

3.8 X 10"^ 

1.6 X 10"^ 
6.4 X 10"" 

4.7 X 10"" 
3.6 X 10"" 
2.4 X 10"" 



b' 


0.46 


1.1 X 10" 


-1 


c' 


1.09 


2.2 X 10" 


-2 


d' 


1.41 


1.2 X 10" 


-2 


e' 


1.94 


3.9 X 10" 


-3 


f 


2.41 


1.8 X 10" 


-3 


g' 


2.85 


5.4 X 10" 


-4 


i' 


3.69 


3.3 X 10" 


-4 


j' 


4.71 


2.3 X 10" 


-4 





Fig. 4. Density contours in the orbital plane of the RLOF type 
flow of Model B (7 = 5/3, Vr = 0.01). The left circle shows 
the mass losing companion star with a radius of 0.25, and the 
right black dot indicates the mass accreting object with a radius 
ofO.015. 



Fig. 6. Comparison of two RLOF type flows. Density contours 
in the orbital plane are shown. Left: Model a (7 = 1.01, Vr = 
0.06); Right: Model D (7 = 5/3, Vr = 0.12). We observe 
spiral shocks which are typical to RLOF. 




0.00 1.31 2.62 3.93 

Fig. 5. Streamlines and Mach number contours in the orbital 
plane of the same model B as in Fig. |2l The Mach number 
contours are represented by the shade, and the numbers in the 
figure show typical Mach numbers. Streamlines are generated 
from points separated by the same distance. We observe that the 
flow rotates counterclockwise around the companion, which is 
a typical feature of RLOF type flows. 



3.3. Intermediate flow 

In the cases E, b and b', we observe that 0.4 < Vr < 0.7. In 
these cases, the flow is not a typical RLOF or a wind accretion 
flow, but rather a very complicated intermediate one. Figure^ 
shows the density contours of model E, and Fig. |8l depicts the 
streamlines of the same model. 




-5. 









1 




-3.33 


.167 


0.00 



Fig. 7. Density contours in the orbital plane of the model E 
(7 = 5/3, Vr — 0.68). We observe very complicated shock 
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0.00 2A 2 4.24 6.36 



Fig. 8. Streamlines and Mach number contours in the orbital 
plane of the same model E as in Fig.0 




Fig. 9. Density contours and stream lines in the orbital plane of 
the model b (7 = 1.01, Vr = 0.46). 



Figure |5] shows the density contours and the streamlines of 
the case b. We find that the flow pattern and the shock structure 
are very complicated in these intermediate flow. 

We will now speculate on the possible cause of the com- 
plicated structure of the intermediate flow. Gas flows from the 
mass-losing star (left sphere) towards the mass accreting object 
(right small sphere). The lower part of the flow experiences a 
shock, and is decelerated to subsonic speed. A small accretion 
disc is formed around the companion, but some part of the gas 
is accelerated to supersonic velocity again. It collides with the 
gas coming from the positive y side and forms a bow shock. 
The orientation of the bow shock is rotated almost 90 degrees 
counterclockwise. 




I 
y 

/ 



Fig. 10. A schematic diagram of an intermediate flow model E. 



3.4. Wind accretion flow 

For the cases with Vr > 0.7, which are all the other models 
than stated above, the gas ejected from the mass-losing star is 
accelerated by thermal pressure. Almost all gas escapes from 
the system and only a small part of gas is accreted by the com- 
panion. We find a typical conical bow shock attached to the ac- 
creting object. Figures^2™d[21show the density distribution, 
the streamlines and the magnitude of the velocity for model g 
(7 = 1.01, Vr ^ 2.86). 




-5,00 -3-33 -1-67 0-00 



Fig. 11. Density contours in the orbital plane of a typical wind 
accretion flow model g (7 = 1.01, Vr — 2.86). We observe a 
conical bow shock attached to the accreting object, which is a 
typical feature of the wind accretion flow. 




0.00 3.62 7.25 10.87 



Fig. 12. Stream lines of the model g as in Fig.fH] Numbers are 
magnitude of normalized velocity. 



3.5. Calculations based on the nested grid 

In order to test the effect of resolution near the mass accreting 
object, we calculate the flow using the nested grid described 
above. Figure [O] shows the density contours for model g'. If 
we compare the result of model g' with that of model g, we 
can see a slight oscillation on the bow shock for model g'. 
Nevertheless, the essential features are unchanged except for 
the mass accretion efficiency. 
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Fig. 13. Density contours of the model g' (7 = 1.01, Vr — 
2.85) calculated on the nested grid. Two blow-ups close to 
the accreting star are shown. Compare the present figure with 
Fig-[n We observe that the bow shock is slightly jaggy com- 
pared with Fig. 10. This is due to the reduced numerical viscos- 
ity. 

3.6. Mass accretion ratio 

We calculate the mass accretion ratio /, which is defined as the 
ratio of the mass accreting rate onto the companion. Mace, to 
the mass loss rate from the mass losing star, Mioss- The mass 
accretion ratio as a function of Vu is shown in Fig. [H] The 
dash-dotted line is an empirical formula / = O.lSx 10^°'^''^ , 
which is valid in the range Vu < 2.5. 

In typical wind case we found that / is less than 1 %, and 
decreases further with increasing Vr although it seems to sat- 
urate (but see below). In the intermediate flow regime, / is 
about a few % to 10%. For RLOF / is greater than 10 % for 
7 — 1.01. If we use the above empirical formula to the limit 
of Vr — 0, we have / = 0.18. Of course, this is not neces- 
sary true and / = 1 is also possible. Although a direct com- 
parison is not possible, our values are in good agreement with 
the results obtained by TBJ96, Walder ( 119971 . Mastrodemos & 
Morris (1998) or Dumm et al. (2000). 

The curve of the calculated accretion ratio based on the 
coarse grid has a kink about Vr = 2.6. This may be due to 
the fact that the size of the accreting cubic box may be com- 
parable to or larger than the accretion radius. In order to see 
the effect, we reduced the size of the inner box. We found that 
the position of the kink is shifted towards larger Vr, and the 
accretion ratio is reduced. This again is in agreement with the 
analysis done by TBJ96. 

For high wind velocity, we expect the Hoyle-Lyttleton the- 
ory of wind accretion (Hoyle & Lyttleton 1939 1, that is the case 
where pressure force is neglected, to be applicable. If it is the 
case, then the mass accretion rate Mace must be proportional 
to the square of the mass accretion radius, 2GM/V^^f - with 
Vinf being the wind velocity at infinity, and so we expect that 
/ ~ Vinf^- Assuming a spherical symmetric wind around the 
mass-losing component, we have 



M^ 



AttA^PaVa, 



(1) 



near the mass-accreting component. For a high wind velocity, 
we may expect that the accretion is determined by Hoyle & 
Lyttleton's formula, i.e.. 



and the accretion radius ra is given by 

_ 2GMi 



(3) 



where pA is the mass density, Va is the wind velocity measured 
on the rotating frame, both at the mass accreting component. 
Ml is the mass of the mass-accreting component. Therefore, 
the accretion ratio is calculated from 



/ = 



Ma, 



wind 



Va 

An 



M- 



Ml + M2 



(4) 



Here, we use Kepler's relation of A^Vl^ — G(Mi + A/2) and 
M2 is the mass of the mass-losing component. For the present 
case of equal mass, the accretion ratio becomes 



/ 



(5) 



Afacc = ■nVaPAVA-, 



(2) 



Note that we checked that the wind velocity does not 
change by a noticeable amount between the Roche lobe sur- 
face of the companion to the position of the primary. Hence, 
we can put Va = Vr in the above equation. 

In Fig. we show the accretion ratio as a function of 
log Vr. The upper line shows the simple Hoyle-Lyttleton for- 
mula as given by Eq. (|5}. The lower line shows another exper- 
imental formula 0.05 Vj^^, which agrees well with the calcu- 
lated results for higher Vr. It appears thus that the calculated 
mass accretion efficiency is about 20 % of the predicted one. 
Such a comparison is however not fair for at least three reasons. 
First, because we take into account pressure, one should not 
use Hoyle-Lyttleton relation but rather the Bondi-Hoyle for- 
mula, and include the sound speed. Second, because our wind 
is thermally driven and we are in a binary system, it is not obvi- 
ous at all to determine which wind velocity should be inserted 
in the equation. In the Hoyle-Lyttleton formaUsm, the velocity 
which is inserted is the upstream wind velocity at infinity. This 
is clearly meaningless in a close binary where the wind may not 
have reached its terminal velocity while close to the compan- 
ion. Third, our wind is a diverging flow from the mass losing 
star, which contradicts the parallel upstream flow assumed in 
the Hoyle-Lyttleton or Bondi-Hoyle treatment. Therefore the 
use of the Hoyle-Lyttleton or Bondi-Hoyle formalism can be 
justified only qualitatively. 

3. 7. Time variation 

We monitored the time variation of the mass accretion ratio. For 
intermediate flows the amplitude of time variation is about 30% 
for 7 = 1.01 and it is about 70% for 7 — 5/3. This trend that 
the amplitude of oscillation is larger for larger 7 is applicable 
to other cases. For larger Vr cases the oscillation is negligible. 
In Fig. [2] the error bars show the amplitude of the oscillation 
of the mass accretion ratio. As was discussed above, the am- 
plitude is larger for the finer grid because of the reduction of 
the numerical viscosity. By decreasing the size of the accreting 
box, we increase the range in Vr showing time variability. We 
Fourier analyzed the time sequence of mass accretion and did 
not find any typical frequency. 
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Fig. 14. Mass accretion ratio, defined by the ratio of the mass 
accretion rate to the mass loss rate, as a function of Vr. The 
horizontal axis is Vr and the vertical one is the mass accretion 
ratio. The dash-dotted line is our empirical formula for the mass 
accretion ratio: 0.18 x 10^"'^^^", while the upper line is the fit 
valid for large Vr. 
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Fig. 15. Log-log plot of the Mass accretion ratio as a function 
of the wind velocity at the position of the Roche lobe. The 
upper line shows a simple Hoyle-Lyttleton formula as given 
by Eq. (|5}- The lower line shows another empirical formula 

/ = 5.0 X 10~^Vj^^. 

4. Conclusion 

1 . The flow pattern is classified by means of the normal speed 
of gas on the critical Roche surface, Vr. For Vr < OAAil, 
we have a RLOF type flow, while for Vr > 0.7 Aft, wind 
accretion flows are realized. In the intermediate parameter 
range, OAAfl < Vr < 0.7 Aft, we observe very complex 
flows. 

2. We construct an empirical formula for the mass accretion 
ratio given by / = 0.18 x iQ-O-TSVn/An^ .^^jj^ f^j. 

3. The calculated mass accretion ratio shows a kink at Vr ~ 
2.6j40. By using the nested grid, we found that this is be- 
cause the assumed size of the mass-accreting box is com- 
parable to or larger than the accretion radius. By reducing 



the size of the mass-accreting box for the nested grid, the 
position of the kink is shifted to larger Vr. 

4. For larger Vr we obtain another empirical formula for the 
mass accretion ratio / — 0.05(VR/yir2)~'''. 

5. For typical wind accretion flow, the results based on the 
coarse grid and those of the nested grid are essentially the 
same except the mass accretion ratio described above. For 
intermediate flow, the computed flow pattern depends on 
the resolution slightly. The amplitude of the oscillation of 
the flow increases with increased resolution because of re- 
duced numerical viscosity. 

Acknowledgements. We thank the anonymous referees for useful 
comments and remarks. Calculations were performed on the Origin 
3800 at the information processing centre of Kobe University and 
on the NEC SX-5 at the Yukawa Institute of Kyoto University. K.O. 
was supported by the Research Fellowships of the Japan Society for 
Promotion of Science for Young Scientists. T.M. was supported by the 
grant in aid for scientific research of the Japan Society of Promotion 
of Science (13640241) and I.H. by (11640226). This work was sup- 
ported by "The 21st Century COE Program of Origin and Evolution of 
Planetary Systems" in Ministry of Education, Culture, Sports, Science 
and Technology (MEXT). 

References 

Anupama, G. C, & Mikolajewska, J. 1999, A&A, 344, 177 
Bennett, P D., Harper, G. M., Brown, A., & Hummel, C. A. 1996, 

ApJ, 471,454 
Biermann, P. 1971, A&A, 10, 205 

Blondin, J. M., Stevens, I. R., & Kallman, T. R. 1991, ApJ, 371, 684 

Blondin, J. M., & Woo, J. W. 1995, ApJ, 445, 889 

Bisikalo, D. V., Boyarchuk, A. A., Chechetkin, V. M., Kuznetsov, O. 

A., & Molteni, D. 1998, MNRAS, 300, 39 
Boroson, B., Kallman, T, Blondin, J. M., & Owen, M. P 2001, ApJ, 

550, 919 

Boffin, H. M. J. & Anzer, U. 1994, A&A, 284, 1026 
Boffin, H. M. J. & Jorissen, A. 1988, A&A, 205, 155 
Boffin, H. M. J. & Zacs, L. 1994, A&A, 291, 811 
Bondi, H., & Hoyle, E 1944, MNRAS, 114, 195 
Carquillat, J.M., Jorissen, A., Udry, S., & Ginestet, N. 1998, A&AS, 
131,49 

Dumm, T, Eolini, D., Nussbaumer, H., Schild, H., Schmutz, W., & 

Walder, R. 2000, A&A, 354, 1014 
Eoglizzo, T, & Ruffert, M. 1997, A&A, 320, 342 
Eoglizzo, T, & Ruffert, M. 1999, A&A, 347, 901 
Eujiwara, H., Makita, M., Nagae, T, & Matsuda, T. 2001, Prog. Theor. 

Phys., 106, 729 

Gawryszczak, A. J., Mikolajewska, J., & Rozyczka, M. 2002, A&A, 
385, 205 

Hachisu, I., & Kato, M. 2001, ApJ, 558, 323-350 
Hoyle, E, & Lyttleton, R. A. 1939, Proc. Camb. Phil. Soc, 35, 405 
Kenyon, S J., Livio, M., Mikolajewska, J., & Tout, C. A. 1993, ApJ, 
407, 81 

Makita, M., Miyawaki, K., & Matsuda, T. 2000, MNRAS, 316, 906 
Mastrodemos, N., & Morris, M. 1998, ApJ, 497, 303 
Matsuda, T, Inoue, M., & Sawada, K. 1987, MNRAS, 226, 785 
Matsuda, T, Ishii, T, Sekino, N., Sawada, K., Shima, E., Livio, M., & 

Anzer, U. 1992, MNRAS, 255, 183 
Matsuda, T, Makita, M., Eujiwara, H., Nagae, T, Haraguchi, K., 

Hayashi, E., & Boffin, H. M. J. 2000, Ap&SS, 274, 259 



Nagae T. et al: Wind accretion in binary stars 



9 



Mikolajewska, J. 1997, in "Physical Processes in Symbiotic 
Binaries and Related Systems," ed. J. Mikolajewska, (Copernicus 
Foundation for Polish Astronomy: Warsaw), p. 3 

Miirset U., & Schmid H. M. 1999, A&AS, 137, 473 

Nussbaumer, H. 1991, Evolution in Astrophysics, Proc. lUE 
Symposium, Toulouse, FR, (ESTEC: Noordwijk) 

Pogorelov, N.V., Ohsugi, Y, & Matsuda, T. 2000, MNRAS, 313, 198 

Ruffert, M. 1996, A&A, 311, 817 

Sawada, K., Hachisu, I. & Matsuda, T. 1984, MNRAS, 206, 673 

Sawada, K., Matsuda, T, & Hachisu, 1. 1986, MNRAS, 221, 679 

Sion, E. M., & Ready, C. 1992, PASP, 104, 87 

Sion, E. M., & Starrfield, S. G. 1994, ApJ, 421, 261 

Sorensen, S.-A., Matsuda, T, & Sakurai, T. 1975, Ap&SS, 33, 465 

Theuns, T. & Jorissen, A. 1993, MNRAS, 265, 946 

Theuns, T, Boffin, H. M. J., & Jorissen, A. 1996, MNRAS, 280, 1264 

Walder, R. 1997, in Accretion phenomena and related outflows, 

Wichramasinghe D.T., Ferrario L., Bicknell G.V. (eds.), lAU 

Colloquium 163, ASP Conf. Ser. No. 121, p. 822 
Walder, R., & Folini, D. 2000, in Thermal and ionization Aspects of 

Flows from Hot Stars, Lamers, H., & Sapar, A. (eds.), ASP Conf. 

Ser. No. 204, p. 331 



